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Dynamics  of  cavitation  clouds  within  a  high-intensity  focused 
ultrasonic  beam 


Yuan  Lu,1  Joseph  Katz,1a)  and  Andrea  Prosperetti1,2 

1 Department  of  Mechanical  Engineering,  Johns  Hopkins  University,  Baltimore,  MD  21218,  USA 
2 Faculty  of  Science  and  Technology,  Impact  Institute  and  J.  M.  Burgers  Center  for  Fluid  Dynamics, 
University  ofTwente,  7500  AE  Enschede,  The  Netherlands 

In  this  experimental  study,  we  generate  a  500  kHz  high-intensity  focused  ultrasonic 
(HIFU)  beam,  with  pressure  amplitude  in  the  focal  zone  of  up  to  1.9  MPa,  in  initially 
quiescent  water.  The  resulting  pressure  field  and  behavior  of  the  cavitation  bubbles  are 
measured  using  high-speed  digital  in-line  holography.  Variations  in  the  water  density  and 
refractive  index  field  are  used  for  determining  the  spatial  distribution  of  the  acoustic 
pressure  nonintrusively.  Several  cavitation  phenomena  occur  within  the  acoustic  partial 
standing  wave  caused  by  the  reflection  of  sound  from  the  wall  of  the  test  chamber.  At  all 
sound  levels,  bubbly  layers  form  in  the  periphery  of  the  focal  zone  in  the  pressure  nodes 
of  the  partial  standing  wave.  At  high  sound  levels,  clouds  of  vapor  microbubbles  are 
generated  and  migrate  in  the  direction  of  the  acoustic  beam.  Both  the  cloud  size  and 
velocity  vary  periodically,  with  the  diameter  peaking  in  the  nodes  and  velocity  in  the 
antinodes.  A  simple  model  involving  linearized  bubble  dynamics,  Bjerknes  forces,  sound 
attenuation  by  the  cloud,  added  mass,  and  drag  is  used  to  predict  the  periodic  velocity  of 
the  bubble  cloud,  as  well  as  qualitatively  explain  the  causes  for  the  variations  in  the  cloud 
size.  The  analysis  shows  that  the  primary  Bjerknes  force  and  drag  dominate  the  cloud 
motion,  and  suggests  that  the  secondary  Bjerknes  force  causes  the  oscillations  in  the  cloud 
size. 

I.  INTRODUCTION 

High-intensity  focused  ultrasound  (HIFU),  along  with  the  associated  cavitation,  is  used  in  a 
variety  of  fields.  The  sonochemistry  involves  the  use  of  acoustic  cavitation  to  initiate  and 
enhance  chemical  reactions.  ’  In  medical  applications,  extracorporeal  shock  wave  lithotripsy 
induced  by  acoustic  cavitation  has  been  investigated  extensively,  both  experimentally  and 
a)  Author  to  whom  correspondence  should  be  addressed.  Electronic  mail:  katz@jhu.edu. 
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numerically,3’4’5’6’7  and  is  used  in  clinical  practice.  HIFU-induced  cavitation  also  plays  an 

o 

important  role  in  delivering  drug  to  target  cells  and  activating  them. 

When  bubbles  are  in  a  non-uniform  acoustic  field,  the  spatial  pressure  gradient  coupled 
with  their  volumetric  pulsation,  produce  a  net  force,  which  is  termed  the  primary  Bjerknes 
force.9  This  acoustically  generated  force  affects  the  bubble  trajectories  along  with  other 
hydrodynamic  forces  such  as  added  mass,  drag,  etc.10'11  During  acoustic  cavitation,  clouds  of 
bubbles  often  form,  either  directly  or  through  fission  of  a  collapsing  bubble,  further 
complicating  the  understanding  and  modeling  of  the  dynamics  involved.  Relative  to  the  single 
bubble  cases,  studies  focusing  on  the  motion  of  cavitation  cloud  in  HIFU  fields  are  relatively 
limited.  Wu  et  al.  irradiated  a  cylindrical  reactor  with  a  27.3  kHz,  low  level  (maximum  sound 
pressure  of  100  kPa)  unfocused  ultrasonic  beam,  and  observed  that  cavitation  clouds  could  be 
generated  in  degassed  water.  The  clouds  traveled  in  the  sound  propagation  direction  toward  the 
water-air  interface  at  a  velocity  of  up  to  1  m/s.  Willard14  employed  an  intense,  2.5  MHz, 
focused  acoustic  beam,  with  pressure  amplitude  of  up  to  7  MPa,  to  study  the  resulting  cavitation. 
He  observed  generation  of  plume-like  bubble  clouds  moving  at  velocity  of  up  to  10  m/s  in  both 
aerated  and  degassed  water.  The  author  suggested  that  the  beam  generated  high-speed  acoustic 
streaming  in  the  liquid,  which  carried  the  clouds.  He  did  not  measure  the  liquid  speed  directly, 
but  inferred  it  from  the  motion  of  microbubbles  around  the  clouds.  The  measured  velocity  did 
not  exceed  2  m/s,  implying  that  other  mechanisms  must  be  involved  in  accelerating  the  clouds. 
The  present  of  neighboring  boundaries  with  the  associated  reflected  waves  further  complicates 
the  dynamics  of  bubble  clouds. 

In  the  present  study,  we  generate  a  HIFU  beam  in  an  otherwise  quiescent  water  container 
with  pressure  amplitude  of  up  to  1.9  MPa,  and  use  high-speed  digital  in-line  holography  to 
measure  the  size,  spatial  distribution,  and  velocity  of  the  bubbles.  Furthermore,  holography 
enables  us  to  record  the  spatial  variations  in  refractive  index  of  water  subjected  to  the  high 
intensity  sound  field,  thus  demonstrating  a  novel  experimental  technique  to  visualize  and 
quantify  the  acoustic  in  the  test  tank.  The  phenomena  we  describe  are  complex.  In  order  to 
better  understand  them,  we  develop  a  model  involving  Bjerknes  forces  and  attenuation  of  the 
sound  field  by  the  cloud.  This  model  is  able  to  predict  the  magnitude  and  spatial  variations  in 
cloud  velocity  reasonably  well. 
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II.  EXPERIMENTAL  SETUP  AND  PROCEDURES 


A.  Facility  and  Instrumentation 

As  illustrated  in  Fig.  1,  the  experiments  have  been  performed  in  a  rectangular  glass  tank 
with  dimensions  of  202x202x101  mm  containing  tap  water.  The  tank  is  open  to  the  air  on  the 
top,  and  the  water  is  initially  quiescent.  The  HIFU  transducer,  manufactured  by  Sonic  Concepts, 
is  positioned  at  the  center  of  the  top  of  the  tank,  with  its  radiating  surface  immersed  in  water. 
This  custom  designed  transducer  operates  at  a  fixed  frequency  /= 500  kHz,  and  has  a  maximum 
input  power  Pm  =400  W.  The  corresponding  wavelength  X  ,  for  a  sound  speed  c0  of  1500  m/s 

in  water  at  20  °C ,  is  3  mm.  The  transducer  has  a  concave  radiating  surface  with  diameter  c/  of 
33  mm,  which  focuses  the  ultrasonic  beam  to  a  zone  located  64  mm  away  from  its  surface.  The 
focal  zone  is  a  cylindrical  volume  with  a  diameter  df  =  5.8  mm  and  a  length  of  42  mm,  as 

provided  by  the  manufacturer  and  confirmed  by  our  measurements  (see  below).  The  500  kHz 
sinusoidal  signal  is  generated  using  a  function  generator,  model  No.  33220A  manufactured  by 
Agilent  Technologies.  The  signal  is  amplified  by  a  200  W  RF  power  amplifier,  model  No. 
1020L,  manufactured  by  Electronics  &  Innovation,  and  monitored  using  a  power  meter,  model 
No.  22 A,  manufactured  by  Sonic  Concepts.  The  signal  is  fed  into  a  matching  network,  also 
manufactured  by  Sonic  Concepts,  and  then  into  the  transducer.  At  high  power  levels,  the 
temperatures  of  the  transducer  itself  and  the  liquid  in  the  focal  zone  increase  rapidly.  To  prevent 
damage  to  the  transducer  and  generation  of  thermally  induced  Rayleigh-Benard  convection  in 
the  liquid,  the  transducer  was  operated  in  a  pulse  mode,  each  pulse  lasting  less  than  2  seconds. 
The  interval  between  two  subsequent  pulses  was  several  minutes  and  long  enough  to  allow  the 
liquid  to  return  to  the  quiescent  condition.  Every  few  pulses,  the  bubbles  accumulating  on  the 
transducer’s  surface  are  removed  to  eliminate  possible  adverse  effects  on  the  ultrasonic  waves. 

A  substantial  fraction  of  the  data  acquisition  about  the  bubble  cloud  characteristics  as  well 
as  the  sound  field  involves  application  of  digital  holography.  Unlike  conventional  photography, 
a  hologram  is  a  record  of  interference  of  coherent  light  scattered  from  objects  with  a  reference 
beam.  Consequently,  it  not  only  contains  information  on  light  intensity  propagating  from 
illuminated  objects,  but  also  the  phase  of  this  light  field.15  Numerical  reconstruction  of  the 
holograms  at  varying  depths  brings  the  objects,  such  as  particles,  bubbles,  etc.,  into  focus, 
enabling  us  to  measure  their  size,  shape  and  spatial  distribution.  ’  ’  ’  ’  ’  For  the  present 


3 


application,  we  employ  digital  in-line  holography  using  the  optical  setup  illustrated  in  Fig.  1. 
Due  to  the  rapid  motion  of  the  bubble  clouds  (velocities  of  the  order  of  m/s),  the  light  source  is 
a  Q-switched  pulsed  523  nm  (green)  laser,  manufactured  by  CrystaLaser,  which  generates  up  to 
0.1  mJ/pulse  at  a  rate  of  up  to  20  kHz.  The  laser  beam  is  collimated  and  illuminates  the  volume 
of  interest,  mostly  in  the  focal  zone  of  the  ultrasonic  beam.  The  holograms  are  recorded  using  a 
high-speed  CMOS  camera,  Photron  FASTCAM-Ultima  APX,  which  can  record  2000  f/s  at  full 
resolution  of  1024x1024  pixels  and  up  to  30,000  f/s  at  a  resolution  of  256x128  pixels.  The 
magnification  varies  from  1:1,  for  which  the  resolution  is  17  //m/pixcl,  to  5:1.  Numerical 
reconstruction  is  performed  using  in-house  developed  software.15 

B.  Visualization  and  measurement  of  the  HIFU  pressure  field 

This  section  describes  the  method  used  to  visualize  the  pressure  field  and  to  characterize  the 
pressure  in  the  focal  zone  of  the  ultrasonic  beam.  The  pressure  wave  in  the  focal  zone  is  nearly 
one-dimensional  and  can  be  expressed  as 

p(x,t)=  p0  +  p'(x,t).  (1) 

where  po  and  p’  are  the  undisturbed  and  acoustic  pressures,  respectively.  The  acoustic  pressure 
can  be  modeled  as  the  sum  of  incident  p'  and  reflected  p'  waves 

r  in  1  re 

=  p'm  //.  =  p{sin(kx-CQt^J  +  CRp{sin(kx  +  CQt  +  gy  (2) 

Here  k  =  2k  /  X  is  the  wave  number,  co  =  2k  f  is  the  angular  frequency,  p{  is  the  amplitude  of 
the  incident  wave,  and  g  is  a  phase  lag.  To  account  for  the  attenuation  of  the  reflected  wave  and 
the  defocusing  of  the  reflected  wave,  we  multiply  its  amplitude  by  a  reflection  coefficient  CR  , 
0  <  CR  <  1 .  The  magnitude  of  p(  can  be  estimated  from  the  input  power  and  the  efficiency  ( 
T]  =  0.85 )  of  the  HIFU  system,  as  provided  by  the  manufacturer.  The  resulting  acoustic  intensity 
at  the  transducer’s  radiating  surface  is  I  r  =  TjP^j^KdrJ ,  and  the  acoustic  pressure  on  the 

surface  is  =  (2IrsZj'~ ,  where  Z  is  the  acoustic  impedance  of  water.  The  manufacturer 
specifies  that  the  intensity  gain  in  the  focal  zone  is  g=2.51,  and  the  incident  pressure  amplitude 
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in  the  focal  zone  is  therefore  pf  =  gp^  .  For  example,  at  P  =  5  W  ,  pf  =  306  kPa ,  while  at 
P  =  110  W  ,  p  =  1.44  MPa  ,  and  at  P  =  200  W  ,  p.  =  1.94  MPa  . 

in  7  r  t  in  7  r  i 

One  can  also  estimate  the  amplitude  of  the  acoustic  field  from  the  changes  to  the  water 
density.  The  substantial  pressure  fluctuations  in  the  focal  zone  cause  detectable  variations  in  the 
density  and  the  refractive  index  n  of  the  water.  When  a  collimated  light  beam  passes  through  the 
spatially  varying  refractive  index  field,  an  initially  collimated  beam  alters  its  direction  of 
propagation.  By  placing  the  hologram  plane  (Fig.  1),  namely  the  plane  being  focused  on  by  the 
imaging  lens,  at  a  certain  distance  H  from  the  center  of  the  acoustic  focal  zone,  the  holograms 
contain  non-uniform  shadowgraphic  patterns  created  by  the  angular  deviations  of  the  laser 
beam.  A  sample  hologram  showing  the  instantaneous  light  intensity  distribution  for  Pm=l  10  W 

is  presented  in  Fig.  2  (a).  The  instantaneous  intensity  fluctuation  of  the  light  field  is  related  to 

22 

the  second  derivative  of  the  refractive  index  normal  to  the  light  axis 


1  z=H  A 


(3) 


where  /o  is  the  undisturbed  intensity  and  I  is  the  intensity  recorded  at  a  distance  z  from  the 
center  of  the  acoustic  beam  (see  Fig.  1).  Assuming  that  the  acoustic  field,  and  thus  the  density 

and  refractive  index,  vary  mainly  in  the  propagation  direction,  V3/?  is  reduced  to  d2n/  dx1  and 

1  \  ^x2)df  ■  For  n  =  n(p),  d2nl  dx2  =  {dnt  dp){d2p  /  dx2) ,  where  dn/dp  is 

22 

a  constant  that  can  be  determined  from  the  Lorenz-Lorentz  relation: 


n-\  1 
n2  +  2  p 


(4) 


where  C  is  a  constant  that  depends  on  the  wavelength  of  the  laser,  properties  of  the  liquid, 
temperature,  and  pressure.  For  a  523  nm  laser  wavelength  at  20  °C  and  atmospheric  pressure, 
C  =  2.1x  10  4  m3/kg  .23  Finally,  using  a  first  order  approximation,  p' =  p'c2  ,  the  acoustic 
pressure  can  be  expressed  in  terms  of  the  light  intensity 


P'U  ,)= _ ^ _ £H 

P  '  ’  ’  (dn/dp)k2dfH  /0 


(5) 
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Given  the  present  experimental  configuration,  the  acoustic  pressure  is  proportional  to  the  light 
intensity  passing  through  the  focal  zone,  and  can  be  estimated  from  the  distribution  of  light 
intensity  (grey  level)  in  the  holograms.  The  bright  bands  in  the  holograms  indicate  locally  high 
pressure  (p'>  0)  since  the  light  converges,  and  thus  brighter  (/'  >0)  in  regions  of  high 
refractive  index  [see  Fig.  2  (a)]. 

In  a  partial  standing  wave,  there  are  no  true  pressure  nodes  since  there  are  no  points  where 
the  rms  pressure  p  is  zero.  However,  for  convenience,  we  still  choose  to  call  the  regions  with 

minimum  p  as  “nodes”  and  those  with  maximum  p  as  “antinodes”.  For  the  acoustic 

1  rms  r  rms 

pressure  field  defined  in  Eq.  (2),  pnas  =  p(  jo.5[j  +  C\  -  2CR  cos (2 Aw  +  g)  j  .  The  reflection 
coefficient  can  be  obtained  from  (7  =  [ pmax  -  //""’)/( pmax  +  pmm) .  Using  Eq.  (5),  p  can  be 

R  \r  rms  1  rms  J /  \r  rms  1  rms  )  ^  1  r  rms 

replaced  with  the  rms  grey  level  (7  )  for  each  pixel,  i.e. 


C„  = 


rms  -*  rms 


/max  _l  1 

rms  i  1 


min 

rms 


(6) 


The  distribution  of  7rms  for  Pjn=l  10  W  is  shown  in  Fig.  2  (b),  for  which  /“=  21.7  ,  7'™"=  13.1 , 
based  on  averaging  along  the  central  line  of  several  antinodes  and  nodes,  and  thus  CR  =  0.25 . 
Using  the  measured  CR  ,  fnm ,  df  provided  by  the  manufacturer,  and  Eq.  (5),  the  estimated 
amplitude  of  the  incident  ultrasonic  wave  is  pf  =  1.32  MPa  .  This  result  differs  by  about  10% 
from  the  previously  mentioned  amplitude  calculated  from  information  provided  by  the 
manufacturer  (p  =  1.44  MPa).  Figure  2  (b)  also  confirms  that  the  width  of  the  focal  zone  is 

consistent  with  the  information  provided  by  the  manufacturer  ( d(  =  5.8  mm  ). 


III.  RESULTS 

A.  Annular  bubble  structures 

24 

At  all  sound  levels,  after  being  generated,  the  bubbles  grow  by  rectified  diffusion.-  Such 
gas  bubbles,  however,  cannot  survive  in  regions  where  the  acoustic  pressure  is  too  high." 
Therefore,  these  gas  bubbles  are  mostly  observed  away  from  the  focal  zone  where,  under  the 
action  of  the  primary  Bjerknes  forces,  they  accumulate  in  layers  separated  by  A  /  2 .  Two 
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examples  are  shown  in  Fig.  3,  the  first  one  showing  the  bubbles  arranged  in  a  vertical  plane,  and 
the  second  one  showing  a  top  view  of  the  circular  structure  surrounding  the  focal  zone.  The 
resonant  radius  of  gas  bubbles  at  the  transducer’s  frequency  and  atmospheric  pressure  is  5.9  pm. 
This  radius  is  smaller  than  that  of  the  bubbles  in  Fig.  3  (a)  implying,  based  on  the  linear  theory 
of  the  primary  Bjerknes  force,  that  the  bubbles  in  this  image  are  attracted  to  the  nodes  of  the 
partial  standing  wave.9  Furthermore,  the  secondary  Bjerknes  forces  among  bubbles  are 
attractive.26  Therefore,  adjacent  bubbles  tend  to  move  toward  each  other,  coalesce,  and  form 
bigger  bubbles.  These  processes  are  evident  in  the  high-speed  movies,  and  the  larger  bubbles  in 
Fig.  3  (a),  which  vary  in  size  between  100-250  pm,  are  the  result  of  the  coalescence  of  several 
small  ones.  The  secondary  Bjerknes  force  also  causes  bubbles  smaller  than  about  20  pm  to 
rapidly  jitter  and  orbit  around  bigger  ones.  When  the  sound  is  turned  off,  these  gas  bubbles  rise 
toward  the  free  surface  due  to  buoyancy. 

B.  Cloud  cavitation 

At  high  sound  levels  (pf  >1.2  MPa ),  cavitation  inception  occurs  near  the  axis  of  the 

acoustic  beam,  not  necessarily  in  the  focal  zone,  and  the  resulting  bubble  clouds  travel  at 
velocities  of  up  to  4  m/s  in  the  direction  of  the  incident  ultrasonic  wave.  This  velocity  is  several 
orders  of  magnitude  larger  than  the  acoustically  induced  streaming,  which  was  measured  by 
tracking  seeded  micro-particles.  Two  examples  of  superpositions  of  high-speed  images 
demonstrating  the  formation  and  migration  of  the  cloud  are  presented  in  Fig.  4.  The  inception  of 
cavitation,  included  in  Fig.  4  (b),  involves  explosive  creation  of  a  nearly  spherical  compact 
volume  of  microbubbles  with  size  of  about  200  pm,  and  rapid  expansion  of  this  cloud  in  the 
horizontal  direction  to  form  a  migrating  layer.  The  inception  occurs  over  a  time  scale  much 
shorter  than  that  resolved  by  the  currently  used  fastest  image  acquisition  rate  of  30,000  f/s,  i.e., 
less  than  33  ps  .  Even  under  these  conditions,  the  velocity  of  the  liquid  surrounding  the 
cavitation  cloud  remains  very  low,  except  for  the  regions  directly  affected  by  the  cloud-induced 
motion,  e.g.  in  the  wake  of  the  cloud.  Thus,  consistent  with  the  observations  by  Willard14,  the 
bubble  cloud  motion  is  not  induced  by  acoustic  streaming.  After  a  few  clouds  are  generated,  the 
surrounding  flow  field  is  sufficiently  disturbed  that  the  bubbly  rings  surrounding  the  focal  zone 
shown  in  Fig.  3  (b)  are  destroyed.  Turning  the  sound  off  makes  the  cloud  disappear  in  less  than 
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33  ps  ,  leaving  no  noticeable  traces.  This  feature  indicates  that  the  bubbles  in  the  cloud  contain 
very  little  non-condensable  gas,  whose  dissolution  would  take  a  longer  time.  These  bubbles, 
therefore,  contain  mostly  water  vapor. 

In  general,  the  cloud  size  increases  with  acoustic  pressure  and  its  dynamics  depends  on  its 
diameter  relative  to  the  sound  wavelength.  When  pf  <  1.68  MPa  ,  the  characteristic  diameter  of 

the  cloud  is  smaller  than  A/3.  In  this  case,  the  cloud  expands  while  decelerating  and  contracts 
while  accelerating  periodically.  The  size  variations  are  evident  in  Fig.  4  (a)  and  (b),  while  the 
changes  in  velocity  are  demonstrated  in  Fig.  5  using  two  examples.  Figures  5  (a)  and  (b)  present 
reconstructed  holographic  images  of  the  clouds,  highlighting  the  changes  in  their  size  and 
location  relative  to  the  acoustic  field,  clearly  showing  that  they  slow  down  near  the  nodes  of  the 
partial  standing  wave.  Corresponding  quantitative  information  on  the  cloud  velocity  and 
diameter  is  provided  in  Fig.  5  (c)  and  (d).  Both  examples  confirm  that  the  cloud  velocity  is 
minimal,  and  its  size  peaks,  near  the  nodes.  On  the  other  hand,  once  formed,  the  thickness  of  the 
clouds  fluctuates,  but  does  not  show  a  clear  correlation  with  the  partial  standing  wave  structure. 
Increasing  the  magnification  of  the  holograms  enables  us  to  resolve  the  spatial  distribution  of 
bubbles  within  the  cloud.  In  the  example  presented  in  Fig.  6,  we  “collapse”  a  series  of  adjacent 
reconstructed  planes  onto  a  single  image.  The  typical  radius  of  the  bubbles  r  is  about  20  pm. 

For  pf  >  1 .68  MPa ,  the  diameter  of  some  of  the  clouds  becomes  comparable  and  even 
larger  than  A/2  [see  Fig.  4  (c)  and  (d)].  In  this  case,  after  the  initial  growth  phase,  the  cloud 
speed  and  diameter  do  not  change  significantly.  As  demonstrated  by  an  unreconstructed 
hologram  in  Fig.  7,  when  a  cloud  with  a  diameter  as  large  as  A  develops  in  the  case  of  high 
sound  level,  the  amplitude  of  the  acoustic  wave  on  the  leeward  side  of  the  cloud  is  much  weaker 
than  that  on  its  windward  side.  This  observation  indicates  that  there  is  little  transmission  of 
sound  through  the  cloud.  Furthermore,  the  amplitude  of  the  wave  on  the  acoustic  windward  side 
appears  to  be  higher  than  that  in  the  surrounding  fluid,  suggesting  a  partial  reflection  of  the 
sound  wave,  which  constructively  interferes  with  the  incident  wave.  These  trends  imply  that,  by 
preventing  sound  transmission,  the  partial  standing  wave  is  eliminated  and,  with  it,  the  periodic 
variations  in  the  cloud’s  size  and  velocity  vanish.  To  confirm  this  statement,  we  covered  the 
bottom  of  the  test  chamber  with  fiber  bristles  in  order  to  attenuate  the  sound  reflection.  Indeed, 
irrespective  of  incident  sound  level,  the  periodic  variations  in  cloud’s  velocity  and  size 


8 


disappeared,  supporting  our  hypothesis  that  both  of  them  are  caused  by  the  partial  standing 
wave.  In  the  following  sections,  we  develop  a  model  that  explains  the  periodic  behavior  of  the 
bubble  clouds. 


IV.  A  MODEL  FOR  THE  AXIAL  MOTION  OF  THE  BUBBLE  CLOUD 


The  cloud  resembles  a  pancake-like  mixture  of  vapor  and  water,  with  volume  of  V  =  7tR2S  , 
where  R  and  5  are  its  radius  and  thickness,  respectively.  However,  some  of  its  features  are  more 
easily  estimated  by  approximating  it  as  a  sphere  with  an  equivalent  radius  of  R  =  (35  /  4 R)1BR  . 

The  velocity  of  the  cloud  is  three  orders  of  magnitude  lower  than  the  speed  of  sound  in  water, 
i.e.  the  associated  time  scales  are  very  different.  Consequently,  it  is  reasonable  to  perform  an 
analysis  of  the  cloud  dynamics,  averaging  over  many  acoustic  cycles,  but  allowing  time  to  vary 
on  the  cloud’s  time  scale.  The  force  balance  on  the  cloud  is 


F  +  F  +  F  +  F 

PB  Alt  D  AM 


0  Vx, 

•  m  5 


(7) 


where  the  overbar  denotes  time  averaging,  dot  denotes  time  derivative  with  respect  to  the  cloud 
time  scale,  FPB  is  the  primary  Bjerknes  force,  FMt  is  the  propulsive  force  due  to  sound 

attenuation  in  the  cloud,  FD  is  the  drag,  FAM  is  the  added  mass, 

pm  =  ( 1  -  (f>)  p  +  (p  pv  ~  ( 1  -  <p)  p  is  the  density  of  the  bubbly  mixture,  cp  is  the  time  averaged 

vapor  volume  fraction,  and  p  and  pv  are  the  time  averaged  densities  of  water  and  vapor, 

respectively.  For  the  present  measurements,  is  about  1%.  For  such  a  volume  fraction,  the 

buoyancy  force  is  two  orders  of  magnitude  smaller  than  the  other  contributions  and, 
consequently,  it  is  neglected.  Note  that  the  buoyancy  force  also  acts  in  the  opposite  direction  to 
the  cloud  motion.  With  expressions  for  all  the  forces,  we  obtain  a  second  order  ordinary 
differential  equation  for  the  cloud  displacement,  and  then  solve  it  numerically. 

The  primary  Bjerknes  force  on  a  cloud  is  (see,  e.g.,  Ref.  9) 


PB 


=  -V'd 


(8) 
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where  the  apostrophe  denotes  the  oscillatory  part  of  the  relevant  quantities.  For  simplicity,  we 
assume  small  oscillations  of  the  cloud  size,  and  estimate  V'  by  using  a  spherical  model,  in 

which  R  is  calculated  from  the  linearized  Rayleigh-Plesset  equation" 


R'  +  2BR'  + Cl2  R'  = 


(9) 


where  [5  and  Q  are  the  damping  coefficient  and  resonant  frequency  of  the  cloud,  respectively, 

and  are  functions  of  x.  Substituting  the  pressure  fluctuations  associated  with  the  partial  standing 
wave  [Eq.  (2)],  we  obtain 


R'  =  — 


Pt 


pR  or 


(10) 
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With  V'  ~  4k R2 R'  ,  the  primary  Bjerknes  force  (8)  becomes 
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The  first  term  in  this  expression  is  associated  with  the  standing  wave,  and  oscillates  with  a 
spatial  period  of  XI 2.  Since  Qr  <  co ,  as  shown  below,  this  contribution  is  always  directed 

toward  the  pressure  nodes.  The  second  term  corresponds  to  the  propagating  part  of  the  wave  due 
to  the  partial  reflection  of  the  incident  wave.  This  part  of  the  force  always  points  in  the  direction 
of  the  incident  wave. 

To  obtain  for  the  cloud,  we  use  the  result  of  a  linear  analysis,  which  shows  that  a 

oo 

spherical  bubble  cloud  has  a  series  of  resonance  frequencies 


3y(i-y)  f rA2 

(j-\!2)2K2{r  , 
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where  j=  1,  2,  3,...,  0)r  is  the  resonant  frequency  of  an  individual  bubble,  and  r  is  its  mean 
radius,  which  is  20  pm  for  the  present  measurements.  Although  this  equation  was  originally 
developed  for  clouds  consisting  of  gas  bubbles,  it  can  readily  be  shown  that  it  is  applicable  for 
the  case  of  vapor  bubbles,  provided  that  co  is  the  natural  frequency  of  a  single  vapor  bubble. 

no 

Since  the  lowest  mode  is  dominant, “  we  will  use  Qr  =  Qj .  Note  that  Qr  is  always  smaller  than 

cor. 

The  damping  coefficient  [5  can  be  estimated  by  relating  the  cloud’s  properties  to  those  of 
the  individual  bubbles  through  the  compressibility  of  the  mixture 

Km=(pmci)1~[(l~v)pcl] 1  >  where  cm  is  the  speed  of  sound  within  the  cloud.  By 

definition,  Km  =  -{l/V^dV  /  dp^j~-{3/ R^dR'J  dpj  and  dR'J  dp'  can  be  obtained  from  the 


solution  to  Eq.  (9),  R'e 


where  i  = 


.  By  taking  the  real 


part  [denoted  with  91  {  }  ]  of  dR'J  dp' ,  we  have 


K 
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(13) 


The  complex  sound  speed  in  a  bubbly  mixture  containing  either  gas  or  vapor  bubbles  is  given 

by29 


3  (p 


—2  2 
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where  b  is  the  damping  coefficient  of  an  individual  bubble.  Since  (l  ~(p)pKm 
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on 

The  expressions  of  co  and  b  of  a  linearly  oscillating  vapor  bubble  are  provided  in  Ref.  : 
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where  3{  }  denotes  the  imaginary  part,  a  is  the  surface  tension,  K  is  the  complex 
compressibility  of  a  vapor  bubble 
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and 


G(w)  =  w4  [ 

Jo 


.exp^-^l  +  i^J 
(w  +  s)5 


ds. 


(20) 


In  the  expressions  above,  J  is  the  ratio  of  the  specific  heats,  cs  is  the  specific  heat  of  the  vapor 
along  the  saturation  line,  L  is  the  latent  heat,  T  is  the  mean  temperature  at  the  vapor-water 
interface,  pb  is  the  mean  pressure  in  the  bubble,  dT  /  dp  ~tJ{  Lpy ) ,  D  and  Z>w  are  the 
thermal  diffusivity  of  vapor  and  water,  respectively,  Kw  is  the  thermal  conductivity  of  water, 
and  A p  is  the  difference  between  the  liquid  static  pressure  and  the  saturation  pressure.  For  the 
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current  case,  2 Z)w/((9r)  =  2.2xl04,  which  is  much  smaller  than  unity.  As  a  result,  B  is 

in  -i  i 

negligible  and  the  calculation  can  be  greatly  simplified.  ’  Using  the  current  parameters, 
(9r  =0.12(9  and  £>  =  0.072(9,  and  as  a  result  Qr  [Eq.  (12)]  is  at  least  an  order  of  magnitude 

smaller  than  (9  .  All  the  terms  needed  for  calculating  [5  [Eq.  (15)]  and  FpB  [Eq.  (11)]  are  now 
available. 

Next,  we  estimate  the  propulsive  force  due  to  sound  attenuation  in  the  cloud,  FA  .  The 

OQ 

attenuation  of  the  acoustic  momentum  flux  (or  the  Reynolds  stress)  in  the  cloud  ’  results  in  a 
net  radiation  force  aimed  in  the  sound  propagation  direction 

FAtt  -  7tR2pAu'2.  (21) 


Here  u'  is  the  acoustically  induced  liquid  velocity  and  A  denotes  the  difference  across  the 
cloud.  For  a  linear  acoustic  wave,  the  velocity  fluctuation  on  the  acoustic  windward  side  of  the 

cloud  is  Mwind  =[Ki-Freexp(-aS)]/(pc0)  ,  and  on  the  leeward  side  is 

w'ee  =  [ //  cxp(-aS)-  p'.c]/(pc()) ,  where  a  is  the  attenuation  coefficient  across  the  cloud.  In  this 

estimate,  we  assume  that  the  cloud  is  small  enough  that  the  local  sound  diffraction  by  it  does  not 
affect  the  overall  strength  of  the  reflected  wave.  Substituting  Eq.  (2)  for  the  pressure,  and  time 
averaging,  we  obtain 


(i-^)(i-c;)Pfw 


2Pco 


(22) 


The  attenuation  coefficient  is  calculated  from  the  complex  sound  speed  in  a  bubbly  mixture 


29 


a  =  -ccB 
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(23) 


Using  Eq.  (14),  and  the  present  conditions,  we  obtain  a  =  8.47xl03  1/m  and 
exp(-2aS)  =  0.034  for  a  characteristic  cloud  thickness  of  200  /dm  ,  i.e.,  the  wave  is 

substantially  attenuated  across  the  cloud. 

The  drag  force  on  the  cloud  is 
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(24) 


Fd  =  ~2Cd  P*R2x\ 

where  CD  is  the  drag  coefficient.  For  a  wide  range  of  Reynolds  numbers,  the  drag  coefficient  of 

a  circular  disk  is  essentially  unity.  Finally,  for  the  added  mass,  we  estimate  it  assuming  a 
spherical  cloud  with  an  equivalent  radius, 

F.\  m  =  =  -  f  PkF%  (25) 


Using  the  derived  expressions  for  the  forces,  and  a  change  of  variables  to  X  =x  +  -^  in 
order  to  eliminate  the  irrelevant  phase  lag  g ,  we  obtain  the  following  nonlinear  second  order 
differential  equation 
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To  solve  this  equation,  we  need  the  dimensions  of  the  cloud  whose  characteristics  are  discussed 
briefly  in  the  following  section,  along  with  its  void  fraction.  To  estimate  the  cloud  motion,  we 
use  the  experimental  observation  that  /?(x)  oscillates  nearly  harmonically  in  space  and,  as  a 


result,  we  estimate  it  as 


r(x)=o.s[(r 


R  )cos2  kX  +  (R 

min  /  \  r 


Based  on  the 


average  dimensions  of  30  clouds,  we  use  R  =350  nm  and  R  =100nm  as  the 

characteristic  maximum  and  minimum  radius,  respectively.  The  thickness  does  not  show  a  clear 
correlation  with  the  partial  standing  wave,  and  seems  to  fluctuate  about  some  mean  value.  For 
the  calculation,  we  use  the  average  thickness  of  the  30  clouds,  S  =  200  /mi .  With  the  initial 
conditions  X(0)  =  A  /  4,  X(0)  =  0 ,  since  cavitation  inception  typically  occurs  in  the  pressure 
antinodes,  Eq.  (26)  is  solved  using  the  Matlab  solver  “ode45”.  Figure  8  (a)  shows  the  “steady- 
state”  solution  for  pf  =  1 .44  MPa ,  and  compares  it  to  the  measured  cloud  velocity  averaged 
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over  30  clouds.  The  model  reproduces  the  correct  magnitude  of  cloud  velocity,  but  there  is  a 
phase  difference.  In  the  case  of  a  big  cloud  (R  =1.5  mm,pf  =  1.89  MPa ),  for  which  CR  =  0 , 
the  model  predicts  a  constant  steady-state  velocity  of  1.36  m/s  (after  the  initial  acceleration), 
which  is  12%  lower  than  the  averaged  observed  value  (1.54  m/s). 

Figure  8  (b)  shows  the  variations  of  calculated  forces  normalized  using  the  predicted  mean 
velocity  and  characteristic  mean  radius.  In  the  vicinity  of  the  nodes  the  total  force  is  negative, 
and  thus  decelerating  the  cloud  there.  The  primary  Bjerknes  force  and  the  drag  contribute  to  this 
deceleration,  while  the  attenuation  force  and  added  mass  are  positive.  Although  the  latter  are 
small,  without  the  attenuation  and  added  mass,  the  cloud  would  not  have  enough  inertia  to 
escape  from  the  nodes,  and  would  be  trapped  there.  Near  the  antinodes,  the  primary  Bjerknes 
force  propels  the  cloud. 


V.  MECHANISMS  AFFECTING  THE  BUBBLE  CLOUD  SIZE 

In  the  previous  section,  the  measured  variations  in  cloud  size  were  used  as  an  input.  In  this 
section,  we  discuss  possible  mechanisms  affecting  the  periodic  contraction  and  expansion  of  the 
smaller  cloud.  As  mentioned  before,  the  radius  of  such  clouds  peaks  at  the  pressure  nodes  and  is 
minimal  at  the  antinodes.  Since  the  cloud  rapidly  disappears  when  the  sound  is  turned  off,  the 
bubbles  much  contain  mostly  vapor  with  gas  diffusion  playing  a  minor  role.  Also,  the  pressure 
oscillations  at  the  antinodes  are  higher  than  those  at  the  nodes,  implying  that  the  larger  cloud  at 
the  node  cannot  be  caused  by  excessive  cavitation  there.  A  possible  mechanism  involves 
interactions  among  bubbles,  such  as  the  secondary  Bjerknes  force,26  namely  the  force  on  a 
bubble  resulting  from  the  pressure  field  generated  by  volume  oscillations  of  all  the  other 
bubbles  in  a  cloud. 

It  is  well  known  that  the  secondary  Bjerknes  force  between  two  equal  bubbles,  according  to 
linear  theory,  is  attractive  both  when  the  bubbles  are  driven  above  and  below  the  resonance 
frequency.  The  theory  gives  this  force  in  the  form 


F  = 

SB 


Anl 


=TV1V2’ 


(27) 


where  v'  and  v'  are  the  oscillating  bubble  volumes,  and  /  is  the  distance  between  the  bubbles. 
To  estimate  this  force,  we  assume  that  the  two  bubbles  have  the  same  equilibrium  radius,  which 
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is  larger  than  the  resonant  size,  and  that  they  are  subjected  to  the  pressure  expressed  in  Eq.  (2). 

27 

Solving  the  linearized  Rayleigh-Plesset  equation"  for  each  bubble  to  estimate  the  volume 
oscillations  gives 


F  = 
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(-2CRcosX  +  C2  +  l). 


(28) 
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As  is  evident,  FSB  is  negative,  i.e.,  it  is  an  attractive  force."  However,  its  magnitude  peaks  at  the 

antinodes,  where  the  pressure  oscillation  peaks,  and  is  minimal  at  the  nodes.  Thus,  the  bubbles 
within  the  cloud  are  more  likely  to  form  dense  clusters  in  the  vicinity  of  the  antinodes.  The 
magnitude  of  the  secondary  Bjerknes  force  can  be  of  the  same  order  as  the  primary  Bjerknes 
force,  especially  when  many  bubbles  are  involved.  Thus,  the  influence  of  FSB  is  not  negligible, 

and  its  variations  along  the  path  of  the  cloud  could  cause  the  changes  to  the  cloud  volume 
observed  during  experiments. 

SUMMARY  AND  CONCLUSIONS 


In  this  study,  we  have  used  high-speed  digital  in-line  holography  to  observe  and  measure 
cavitation  phenomena  within  a  high  intensity  focused  ultrasonic  beam  in  otherwise  quiescent 
water.  The  amplitude  of  pressure  oscillations  is  sufficiently  high  to  cause  corresponding 
changes  to  the  water  density  and,  consequently,  to  its  refractive  index.  These  changes  are 
recorded  in  holograms  and  used  to  estimate  the  amplitude  of  pressure  oscillations.  The  partially 
attenuated  wave  reflected  by  the  container  walls  interference  with  the  primary  wave  giving  rise 
to  a  partial  standing  wave,  in  which  the  pressure  field  consists  of  both  traveling  and  standing 
components  with  clearly  defined  pressure  nodal  and  antinodal  regions.  Clouds  of  cavitation 
bubbles  typically  appear  first  near  the  antinodes,  and  then  travel  at  speeds  of  several  m/s  in  the 
same  direction  as  the  primary  sound  wave,  consistent  with  prior  observations.13'14  The  typical 
cloud  size  is  a  fraction  of  a  millimeter,  much  smaller  than  the  ultrasonic  beam  size,  and  it 
contains  hundreds  of  microbubbles  with  typical  radius  of  20  pm.  Since  the  bubbles  disappear  in 
less  then  our  fastest  frame  rate  of  33  /ils,  we  conclude  that  they  contain  mostly  vapor.  Therefore, 
the  diffusion  of  non-condensable  gases  dissolved  in  the  liquid  does  not  play  a  significant  role. 
The  size  and  speed  of  the  smaller  clouds  vary  periodically,  with  the  speed  peaking  at  the  nodes, 
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and  the  size  at  the  antinodes.  When  the  cloud  becomes  large  enough  to  block  a  significant 
fraction  of  the  acoustic  beam,  these  variations  disappear.  Under  such  conditions,  after  an  initial 
growth,  the  cloud  size  and  1-2  m/s  propagation  velocity  remain  approximately  constant. 

The  cloud  motion  cannot  be  a  result  of  acoustic  streaming  since  the  bubbles  are  transported 
at  speeds  that  are  much  higher  than  that  of  the  liquid  surrounding  the  cloud.  To  explain  the 
observed  phenomena,  we  have  used  a  model  that  accounts  for  the  primary  Bjerknes  force,  sound 
attenuation,  added  mass  and  drag.  The  volumetric  oscillation  of  the  cloud,  required  for 
calculating  the  primary  Bjerknes  force,  as  well  as  the  sound  attenuation  through  the  cloud,  are 
estimated  using  linearized  theories.  With  the  cloud  size  as  an  input,  the  model  predicts  fairly 
well  the  velocity  magnitude  and  its  variations  along  the  partial  standing  wave.  The  dominant 
forces  in  the  cloud  dynamics  are  the  primary  Bjerknes  force  and  the  drag  force.  The  model  also 
predicts  the  cloud  motion  when  the  sound  attenuation  is  high  enough  to  eliminate  the  reflected 
wave. 

The  oscillations  in  cloud  size,  which  peaks  in  the  nodal  regions,  cannot  be  caused  by  spatial 
variations  in  cavitation,  since  the  pressure  oscillations  peak  at  the  antinodes.  To  explain  them,  at 
least  qualitatively,  we  consider  the  effect  of  the  secondary  Bjerknes  force,  which  involves  the 
interaction  among  bubbles.  The  force  calculated  using  linear  bubble  dynamics  indicates 
attraction  among  the  bubbles,  which  increases  with  the  magnitude  of  pressure  oscillations  and 
can  be  of  the  same  order  as  the  primary  Bjerknes  force.  This  attraction  would  reduce  the  cloud 
size  around  the  antinodes  to  a  greater  extent  than  in  the  vicinity  of  the  nodes,  possibly 
explaining  the  observed  trends. 
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FIGURES  AND  CAPTIONS 
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Fig.  1.  (Color  online)  Setups  for  the  acoustic  cavitation  and  high-speed  digital  in-line  holography  measurements. 
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Fig.  2.  Visualization  and  quantification  of  the  focal  zone  of  the  ultrasonic  wave;  pf=\AA  MPa.  (a)  A  sample 
hologram  showing  the  instantaneous  acoustic  wave.  The  bright  bands  correspond  to  high  pressure,  (b)  The 
distribution  of  the  rms  grey  levels,  showing  the  partial  standing  wave  structure. 
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Fig.  3.  (a)  Annular  bubble  layers  in  the  pressure  nodes  at  low  sound  intensity  (pf=306  kPa)  recorded  using  white 
light  illumination,  (b)  A  top-view  of  an  early  test  chamber  showing  that  the  bubbles  accumulate  in  the  periphery  of 
the  focal  zone  (denoted  with  dashed  circle). 
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Fig.  4.  Superpositions  of  time  series  of  sample  bubble  clouds  recorded  using  white  light  at  (a-b)  high  sound  level 
{p t=1.44  MPa),  where  the  speed  and  diameter  of  the  cloud  are  periodic,  with  (a)  showing  a  few  samples  and  (b) 
containing  the  entire  series;  (c-d)  very  high  sound  level  (pt=1.89  MPa),  where  the  cloud  becomes  large,  travels  at 
almost  a  constant  velocity,  and  maintains  constant  diameter  after  the  initial  growth  phase.  Here  (c)  shows  a  few 
samples  and  (d)  the  entire  series. 
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Fig.  5.  (a-b)  Time  series  of  reconstructed  holograms  (Af=167  /us)  showing  the  axial  locations  of  two  sample  bubble 
clouds  in  a  motion  cycle,  and  the  corresponding  location  of  the  pressure  nodes  (thick  solid-lines)  determined  from 
the  distribution  of  rms  pressure  levels;  pf=  1.44  MPa.  (c-d)  Corresponding  velocity  and  diameter  vs.  the  cloud 
location  in  a  motion  cycle. 


25 


980  pm 


Fig  6.  A  close-up  view  of  the  inner  structure  of  a  bubble  cloud;  pf=  1.44  MPa.  This  picture  is  generated  by 
collapsing  a  series  of  reconstructed  holograms  from  different  depths  onto  a  single  plane. 
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Fig.  7.  A  sample  original  instantaneous  hologram  showing  the  reflection  of  the  incident  acoustic  wave  by  a  large 
bubble  cloud;  pf=1.89  MPa.  The  wave  on  the  leeward  side  of  the  cloud  becomes  nearly  invisible.  The  bubbles  that 
are  out  of  focus  can  be  seen  in  the  background. 
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Fig.  8.  (a)  A  comparison  between  the  measured  and  predicted  [Eq.  (26)]  bubble  cloud  velocity.  The  error  bars 
indicate  the  standard  deviation  of  the  measured  values,  (b)  The  predicted  total  force  and  its  components;  p(=  1.44 
MPa. 
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